Superconductivity in a magnetically ordered background 
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Borocarbide compounds with the formula RM2B2C show 
interesting superconducting and magnetic properties and the 
coexistence of the two phenomena. BCS theory is extended 
to systems with underlying commensurate magnetic order. 
In the case of helical phases the technique may be extended 
to any Q-vector and there exists a well defined limit for in- 
commensurate values. The way magnetic order influences su- 
perconductivity depends crucially on the details of both the 
magnetic structure and the electron bands, but some quali- 
tative criteria may be given. As an example we give a brief 
analysis of the compound HoNi2B2C. 

The borocarbides are a class of compounds with for- 
mula RM2B2C (R = Y, La, or a rare earth). They at- 
tract attention because of their interesting superconduct- 
ing and magnetic properties and their mutual interaction. 
Compounds with Y, Lu, Tm, Er, Ho and Dy have phonon 
mediated superconductivity at relatively high tempera- 
tures Q. On the magnetic side Gd, Tb, Dy, Ho, Er and 
Tm show long-range magnetic order with a variety of 
different structures ||. Crystalline electric field effects 
determine the easy direction of the spins in the localised 
/-states of the rare earth ions, while the RKKY inter- 
action among them determines the ordering Q vector. 
Commensurate and incommensurate magnetic structures 
coexist with superconductivity. The case of HoNi 2 B 2 C 
has been the most spectacular, where a c-axis incom- 
mensurate helix, an a-axis incommensurate phase and a 
c-axis antiferromagnet appear in the same temperature 
region in which the superconducting upper critical field 
shows a pronounced depression 

For the description of superconductivity in the coex- 
istence region we use an extended BCS theory. The 
Hamiltonian of band electrons interacting with a given 
exchange field is: 

H M = £k C k* CkCT ~ 12 h 1 ' <T<T<T ' C k+qa C ka' , (1) 
ko" kqerer' 

where h q is the Fourier transform of the internal effective 
magnetic field (h q = [ib 9 -f (q) S(q)) and a Ga is a vector 
with the Pauli matrices as the components. Conventional 
notation is used for the other symbols || . This has to be 
added to the usual phonon Hamiltonian: 

HsC = y^tJq ajflq + 12 5k ' k C kV C k<r (ttk'-k + Ok_k') 
q k'ker 

(2) 

For a given exchange field h q commensurate with the lat- 
tice the eigenstates problem for Hm is reduced to solving 



a nxn matrix, where n is the number of rare earth atoms 
in the magnetic unit cell. Using the extended unit cell 
representation, the total Hamiltonian of the system is 
given by: 

H = ^ ? ki/ c^cW +^]w q a+a q + 

\tu q 

+ 12 9k<kc£>v'Ckv (ah/_k + a£_ k ,) (3) 

k'kz^V 

The main differences with respect to the non-magnetic 
case are: 

(1) the magnetic eigenstates are labeled by the index v 
that does not refer to a definite spin state. Furthermore 
energy spin degeneracy is in general lifted and magnetic 
Cooper pairs formed with these states will not be in spin 
singlet state 

(2) the band structure acquires gaps at magnetic Bragg 
planes and possibly at other characteristic surfaces. 
Whenever the Fermi surface intersects one of those planes 
or surfaces it is force to become orthogonal to it. How- 
ever, apart from pathological situations, the properties 
of the electrons are affected only to order — 

(3) the scattering amplitude between magnetic states 
5k'k acquires further k dependence linked to the under- 
lying exchange field. Typical magnetic energy bands are 
shown in the figure . 

Introducing the gap function matrix: 

A yV = (cWc-k,) (4) 

we are able to construct the mean-field theory of super- 
conductivity in a way quite close to the standard one. 
Some qualitative remarks are now possible: in mag- 
netic states with a net magnetisation (i.e. ferromag- 
nets) v = a, spin degeneracy is lifted and coexistence 
is possible only for very small values of the magnetisa- 
tion. On the other hand, if the net magnetisation is zero 
and the magnetic structure is collinear to a given vec- 
tor n (S(x) = 5(x) • n), the spin along that direction is 
conserved. In these cases (i.e., antiferromagnets |jj and 
longitudinally modulated phases) the magnetic theory re- 
duces to the usual one but with modified band structures 
and interaction with phonons. 

At a qualitative level it may be noticed that: 

(1) if magnetic order does not destroy large pieces of the 
Fermi surface, as in the case of nesting, electronic prop- 
erties are only slightly modified and the opening of the 
magnetic gap does not imply strong suppression of su- 
perconductivity 

(2) time-reversal symmetry of the system is not a neces- 
sary condition for the existence of superconducting pairs, 
nor is time-reversal plus a translation symmetry 
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Now we consider the case of helical magnetic order. For 
a helix with wave vector Q and amplitude h, the mag- 
netic eigenstates may be computed analytically through 
a Bogoliubov transformation: 



e k±Q + £k £k — £k±Q 



4/l 2 
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and the Bogoliubov coefficients are: 



(ek - £k±Q) 2 
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(fk - £k+Q) 2 
(£k - ek+q) 2 + 4fc 2 



« 2 (k) = 1 - u 2 (k) 



Cooper pairs formed with these electronic states are in a 
mixed state of singlet and triplet. Using the Einstein ap- 
proximation for the phonon dispersion law (u>q) and the 
structureless electron-phonon interaction (<?k'k = g) in 
the weak coupling limit it is possible to obtain the stan- 
dard BCS equation for the gap function with a modified 
definition of the interaction parameter (see Ref. for 
more details): 



A = 



w 



dS 



(u 2 (k)-v 2 (k)Y 
|V k e k | 
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In order to give a quantitative estimate for the new effec- 
tive electron-phonon coupling constant we assume cylin- 
drical symmetry around the z-axis for the energy bands 
and linearise them at the magnetic gap. In this way 
band effects enter through the two components of the 
Fermi velocity (fivp = V k e k ): the one orthogonal (v±) 
and the one parallel (f|| ) to the z-axis. The difference in 
the interaction parameter is given by: 



AA = A - A 
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where k r ~ kp is the radius of the intersection between 
the Fermi surface and the magnetic Bragg plane. The 
relative reduction of the superconducting parameter is 
then: 



AA 



(%Vpk r ) h 



(hv±kp) (hvukp) 



(5) 



Eq. (g) has the following qualitative features: 

(1) assuming v± and Vu of order vf the variation of the 
interaction parameter is proportional to j- and A is re- 
duced only by a few percent 

(2) however RKKY interaction maxima are usually con- 
nected to pieces of the Fermi surface orthogonal to the Q 
vector close to the Bragg planes. If the Q vector of the 
magnetic order is determined by the RKKY interaction 
(as it happens generally for incommensurate structures), 



then v± <C vp and the suppression of superconductiv- 
ity may be large. These qualitative features may give 
an interpretation of the behaviour of superconductivity 
in HoNi2B2C. From anisotropic magnetisation measure- 
ment a broad maximum in the RKKY function is inferred 
close to q = (0, 0, n) || . This maximum is related to both 
the magnetic c-axis incommensurate helix and to the o 
axis antiferromagnet. Because of its flatness it is prob- 
ably due to the overall structure of the electronic bands 
and not to sharp nesting features, thus v± ~ vp. On the 
other hand the a-axis modulation appears to be related 
to nesting features of the Fermi surface (both experimen- 
tally and through band calculations ||). We therefore 
suggest that the a-axis modulation has a stronger de- 
structive character for superconductivity than the c-axis 
modulation Jl0| . 

We conclude that an extended BCS theory may be con- 
structed in the general case of commensurate magnetic 
structures. In the special case of helical magnetic order 
the extended theory is also valid in the incommensurate 
limit. A possible qualitative explanation for the anoma- 
lies in the superconductivity of HoNi 2 B 2 C has been sug- 
gested. 
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FIG. 1. Effect of the magnetic order on two model band 
structures. The magnetic vector used are q = ir in the middle 
row, corresponding to antiferromagnetic order, and q = 4^ in 
the lower row, corresponding to structure with periodicity of 
three unit cells. Note that magnetic gaps here correspond to 
~ 0.1 in order to emphasise them, physical values for the 
gaps are 10 to 100 times smaller. 



